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Abstract. Lie symmetry analysis is applied to study the nonlinear rotating shallow 
water equations. The 9-dimensional Lie algebra of point symmetries admitted by the 
model is found. It is shown that the rotating shallow water equations are related with 
the classical shallow water model with the change of variables. The derived symmetries 
are used to generate new exact solutions of the rotating shallow equations. In 
particular, a new class of time-periodic solutions with quasi-closed particle trajectories 
, is constructed and studied. The symmetry reduction method is also used to obtain 

some invariant solutions of the model. Examples of these solutions are presented with 
a brief physical interpretation. 

7— I ■ 
> . 

O PACS numbers: 92.05.Bc, 92.10.Ei 

ON 

oo 

o 

oo 

o 



AMS classification scheme numbers: 35C05, 76B15 



1. Introduction 



The article focuses on the Lie symmetry analysis of the rotating shallow water (RSW) 
model. In Cartesian frame of reference the RSW equations have the following form 

u t + uu x + vu y - fv + gh x = 0, 

v t + uv x + vv y + fu + ghy = 0, (1.1) 

ht + {uh) x + (vh) y = 0. 

Here (it, v) is the fluid velocity, h is the free surface height over the flat bottom, / is the 
constant Coriolis parameter and g is the constant gravity acceleration. The RSW model 
(11. ip arises from the three-dimensional rotating incompressible Euler equations (see p], 
pp. 59-63) under the following assumptions: a) the scale H of the vertical motion is 
much less than the scale L Q of the horizontal motion, so that H /L << 1; b) the fluid 
density is constant (p = const); c) the external force is due to gravity, and pressure 
obeys the hydrostatic approximation p = pg(h — z) + p (p is the constant pressure on 
the free surface); d) the axis of rotation of the fluid coincides with the vertical z-axis. 

This nonlinear system of partial differential equations (PDEs) is widely-used 
approximation for atmospheric and oceanic motions in the midlatitudes with relatively 
large length and time scales. The model applies to phenomena that do not depend 
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substantially on temporal changes of the density stratification. In spite of its simplicity, 
it contains all essential ingredients of atmosphere and ocean dynamics at the synoptic 
scale. The main known mathematical results and physical applications concerning of 
the RSW model (11.11) as well as its linear theory, quasi-geostrophic limit and etc are 
presented in PQ [3]. 

We recall here that the potential vorticity, defined by Q = (v x — u y + is 
conserved on particles: Q t + uQ x + vQ y — 0. Moreover, there are two important for 
applications non-dimensional parameters: Rossby number Ro = Uo/(fLo) and Froude 
number Fr = Uq / y/gHo (here Uq is the typical velocity scale for the flow). In non- 
dimensional form of the model Ro -1 and Fr~ 2 stand for / and g, correspondingly. In 
this study these parameters do not play a vital part. In fact, as result of obvious scaling 
transformations the system ( 11. II) can be reduced to the equivalent one with / = 1 and 
g = 1 (in case / ^ and g ^ 0). 

One of the methods for studying of PDEs is group analysis. This analytical 
approach based on symmetries of differential equations was originally introduced by 
Sophus Lie at the end-nineteenth century and further developed by Ovsynnikov [5], 
Olver [6] and others. For each system of PDEs there is a symmetry group, that acts on 
the space of its independent and dependent variables, leaving the form of the system 
unchanged. The classical Lie symmetry analysis allows one to construct and study 
special types of analytical solutions of nonlinear PDEs in terms of solutions of lower- 
dimension equations. For modern description of the theory see, for example, [5]-[7]. 
Many applications of group analysis to PDEs are collected in [8]. The classical Lie 
method is an algorithmic procedure for which many symbolic manipulation programs 
were designed (SlUOj- This software became imperative in finding symmetries associated 
with large systems of PDEs. 

Exact solutions of nonlinear systems descriptive of fluid motions with moving 
boundaries are uncommon even in the shallow water approximation. Freeman [UJ and 
Sachdev [12] have obtained solutions for gravity wave motions under the assumption that 
the flow is a simple wave. It is interesting to note that the similarity solution obtained 
in [11] may be systematically generated by group methods [13]. In the latter work, 
symmetry analysis was often used to obtain analytical solutions to the nonlinear water- 
wave problems. Symmetries and exact solutions of a nonlinear system that models the 
finite motion of a rotating shallow liquid contained in a rigid basin were studied in [H] 
(see also [15]). At present there known many examples of exact (invariant and partial 
invariant) solutions for nonlinear models of fluid dynamics [13] [27] and others. Amount 
of analytical solutions to the Euler equations govern the compressible inviscid flow were 
obtained within the framework of Ovsynnikov's "submodels programme" [T6] [T7]. 

In this paper we investigate Lie point symmetries and classes of exact solutions 
of the RSW equations. In Section 2, we determine the 9-dimension Lie algebra of 
infinitesimal transformations admitted by the RSW model. We also prove that the 
derived Lie algebra of symmetries is isomorphic to the Lie algebra of infinitesimal 
transformations admitted by (2+l)-dimensional shallow water (SW) equations. This 



Symmetries and exact solutions of the rotating shallow water equations 



3 



allows one to use its known optimal system of subalgebras. Moreover, we find the change 
of variables which transforms the RSW model to the SW equations. This is one of the 
principal results of the paper. In Section 3, we derive and analyse finite transformations 
corresponding to the nontrivial symmetries of the RSW model. On basis of these 
transformations, in Section 4, we construct new time-periodic exact solutions of the 
RSW equations. These solutions may be interpreted as pulsation of liquid volume under 
the influence of gravity and Coriolis forces. In Section 5, using some two-dimensional 
parameterized classes of the optimal system of subalgebras we reduce the RSW model 
to the ordinary differential equations (ODEs) and integrate them. In particular, we 
construct and study exact solutions describing rotational-symmetric subcritical and 
supercritical regimes of the stationary flow as well as various non-stationary regimes 
of spreading and collapse of a liquid ring. 



2. Symmetries of the RSW equations and reduction to the SW model 

2.1. Lie point symmetries of the RSW equations 

On the basis of the group analysis of the differential equations we shall find the 
symmetry group of infinitesimal transformations associated with the RSW equations 
( II. ip . Following [5], let us define infinitesimal generator X and its first prolongation Y 

X = f(x, u)^ + rf(x, u)oU Y = X + Qd^ = 1, 3). 

i 

Here we use the following notations for convenience 

x = (ar, x 2 , x 3 ) = (t, x,y), u = (u , u 2 , u 3 ) = (u, v, h) 



dx v dx l l dui 

To determine the group of symmetries admitted by the system of equations (11.11) . we 
act on it with the first prolongation of generator X and go to the set of the system 
solutions. Consequently, we get a system of the determining equations for the unknown 
functions £ J (x, u) and ^ (x, u), that allows splitting in variables u\. Let us omit the 
bulky intermediate calculations and present the final result of the symmetry group 
determination. 

Theorem 1 The symmetry group associated with the RSW equations (GUP is generated 
by the following vector fields: 

X\ = d x , X2 = d y , Xj = dt, 

X 3 = cos(ft)d x - sin(ft)d y - f sin( ft) d u - fcos(ft)d v , 

X 4 = sm(ft)d x + cos(ft)d y + fcos(ft)d u -fsin(ft)d v , (2.1) 

X5 = -yd x + xd y - vd u + ud v , 

X 6 = xd x + yd y + ud u + vd v + 2hd h , 
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X 8 = cos(ft)d t - -(xsm(ft) -ycos(ft)^jd x - -(xcos(ft) +yean(ft))d y + 
+|((« " fy) Mft) + (v- fx) cos(/t))9„ - 

"{((«+ /V) cos (^) - ( v + / x ) Mft))d v + fh sm( ft)d h , (2.2) 

X 9 = sm(ft)d t + cos(/t) + y sin(/*)) & - ^ (a; sin(/t) - y cos(ft)j d y - 

~ 2 _ ^ cos ^ - ( u - Z^) sin(/*)^9„ - 

+ /y) sm(/t) + (v + fx) cos(/t))^ - fhcos(ft)d h . (2.3) 

Thus, the RSW equations ( 11. ip are invariant under translations in t, x and ?/ (X7, 
Xi and X 2 , correspondingly), helical rotations with respect to t (X 4 and X 5 ), rotation 
(X 5 ), scaling symmetry (X 6 ) and two complicated vector fields (X 8 and X 9 ) which will 
be discussed below. The infinitesimal symmetries fl2.ip - fl2.3l) form Lie algebra Lg. 

Symmetry properties of the RSW equations in Lagrangian coordinates were studied 
in [28]. Transition to Lagrangian variables is a nonlocal change of coordinate. Therefore 
Lie groups admitted by the model in Eulerian and Lagrangian variables do not 
completely coincide. In particular, there are no nontrivial symmetries defined by 
generators Xg and Xg in Lagrangian variables. 

As mentioned above, in [HJ [15] group analysis was applied to the model of rotating 
shallow liquid in a rigid basin with the following geometry 

z = Z(x, y) = Ax 2 + By 2 , A > 0, B > 0. 

It was shown that in the elliptic paraboloid case, A / B, the symmetry algebra is 
spanned by six infinitesimal generators. In the circular paraboloid case, corresponding 
to A = B, the symmetry algebra is larger — namely, 9-dimensional. Note that these 
symmetries reduce to (I2.1l) - fl2.3p for A = B = 0. 



2.2. Properties of Lie algebra L 9 

Let E be a given system of differential equations admitting a symmetry group G. The 
basic property of solutions for E is that any solution of the system E carried over by any 
transformation of the group G to a certain solution of the same system E. Therefore, two 
solutions of the system E are said to be essentially different with respect to G if they can 
not transformed to each other by any transformation of the group G. Thus, it is useful 
to enumerate subgroups of the group G, which lead to essentially different solutions. By 
virtue of the fact that there are certain correspondences between the symmetry groups 
and Lie algebras, it is enough to construct an optimal system of subalgebras QL (see 
[5], [291 E0]), i-e. the minimal set of subalgebras of L, which exhaust all the essentially 
different invariant and partially invariant solutions to equations E. 

In the present case there is no need to construct the optimal system of the 
subalgebras for Lie algebra Lg since we can use the results derived from the symmetry 
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analysis of the gas dynamic equations. In the case of zero Coriolis parameter (/ = 0) the 
model under consideration reduces to the SW equation (or two-dimensional polytropic 
gas dynamic equations with polytropic exponent 7 = 2) 

u t + uu x + vuy + gh x = 0, 

v t + uv x + vv y + gh y = 0, (2.4) 

ht + (uh) x + {vh) y = 0. 

The Lie algebra Lg of infinitesimal symmetries of the SW equations is spanned by the 
following vector fields [HI [8] 

Z\ = d x , Z 2 = d y , Z 3 = td x + d u , Z4 = td y + d v , 

Z 5 = -yd x + xd y - vd u + ud v , 

Z 6 = xd x + yd y + ud u + vd v + 2hd h , Z 7 = d t , 

Z 8 = t 2 d t + txd x + tydy + (x — tu)d u + (y — tv)d v — 2thd h , 

Z 9 = 2td t + xd x + ydy — ud u — vd v — 2hd h . 

For this Lie algebra Lg, the optimal system of subalgebras QLg, which includes 
179 parameterized classes is presented in [18] . By virtue of the following theorem, 
parameterized classes of this optimal system of subalgebras can be applied for obtaining 
invariant or partially invariant solutions of the RSW equations (11. lft . 

Theorem 2 Lie algebras Lg(Zi, Zg) and Lg(Xi, Xg) are isomorphic. 

Proof. Let us redefine the basis of generators of the Lie algebra Lg so that 
Yi = X 2 — -X4, Y 2 = X 3 — Xi, Y 3 = X 1 + X 3 , Y A = X 2 + X4, 



Ys — X 5 , Y 6 — X 6 , Y 7 — — ( X.j — X 5 — X t 



Y& — — (^Xj — -^Xs + X 8 j , Yg — ——Xg. 

The commutation relations for the Lie algebra Lg in this basis are displayed in table 
CD The basis is a canonical one chosen so that the Levi decomposition is immediately 
apparent from the commutation table. Thus Lg = R(B N, where R = {Yi, Y 6 } and 
N = {Yj,Ys,Yg} constitute the maximal solvable ideal (radical) and the simple Lie 
algebra sl(2). The nilradical of Lg is abelian and is generated by {Yi, Y 2 , Y 3 , Y 4 }. 

It is easy to verify that the commutation relations for the Lie algebra Lg(Zi, Zg) 
are defined by the same table of commutators (where symbols Y/c should be replaced 
by Zk). The table of commutators for Lie algebras L 9 (Zi, Zg) and Lg(Yi, Yg) are 
coincide completely. Consequently, Lie algebras Lg and Lg are isomorphic. □ 

Thus, to construct invariant or partially invariant solutions to the RSW equations 
( 11.11) we can use generators Yk defined by formulas ( I2.1I) - (12.3I) and (12.51) as well as the 
parameterized classes of the optimal system of the subalgebras introduced in 1181 . 
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Table 1. Table of commutators for Lie algebra Lg(¥x, ...,Yg). 
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-Y 4 
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Y 6 


-Yt 


-Y 2 


-Y 3 


-Y A 

















Y 7 








Y x 


Y 2 











Y 9 


2Y 7 


Y 8 


-3^3 


-Y 4 














—Yg 





-2Y 8 


Y 9 


-Y t 


-Y 2 


Y 3 


Y 4 








-2Y 7 


2Y 8 






^.<?. Reduction to the SW model 

In view of the previous theorem arises important question. Is there change of variables 
such that each generator {k = 1, ...,9) reduces to generator Z^! In the case of 
existence one it is necessary to verify if the RSW equations (11.11) reducible to the SW 
equations (12 .4p under this transformation. The answer on the question is given by the 
following theorem, which can be easy proved by straightforward calculations. 

Theorem 3 The models U.l\) and \2.J$ as well as their symmetries are related by the 
following transformation 

< / = ~oot(-^), ^ = -l(xcot{^)-v), y' = -l(x + ycot(^)), 

u' = - l -{u sm(/t) -v(l- cos(/t)) - fx) , (2-6) 

v' = ~\ {u{l - cos(/t)) + v sm(ft) - fy), h' = ± (l - cos(/t)) . 

In fact, let us suppose that a set of functions u(t,x,y), u(t,x,y), h(t,x,y) satisfies 
to the RSW equations (II. ip . Then the immediate calculations show that the following 
set of functions u'(t', x', y'), v'(t', x', y'), h'(t', x', y'), given by formulae (I2.6p . satisfies to 
the SW equations (I2.4p . Thus, the models (ll.ip and (12 .4p are related by the change of 
variables (I2.6p . 

It is easy to clear up that each generator Y& in the variables with prime (12.61) takes 
form Z k . More precisely, Yj, = Z k (k = 1, 2, 5, 6, 9), Yj, = fZ k (k = 3, 4, 8), Y 7 ' = f~ x Z 7 . 
Using scaling it is possible to fix formulae (I2.6P (or reduce (11.11) to the equivalent one 
with / = 1) in order to Y" fc ' = Z k for all k — 1, 9 without coefficient /. However, this 
form is more convenient as it conserve physical dimension of the variables. 

Note that the mapping (12. 6p transforms any solution of the SW equations (12.41) 
defined for all time moments (— oo < t < oo) to the solution of the RSW equations (11.11) 
defined only for bounded time interval (0 < t < 2n/ f) or, may be, to the stationary 
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rft\ fr rft\\ , ^ 

u -v cot^—J--{x-ycot^—JJ, h = Y 



solution. This change of variables allows one to construct solutions of the RSW model 
using known solutions of the SW equations and vice versa. 

For example, let us consider constant solution of the SW equations ( 12 Ah . u = Uq, 
v = vq, h = ho- Obviously, in this case trajectories of fluid particles are lines. According 
to ( 12.61) we obtain the following solution of the RSW model (11. ip 

u = u cot - v + - (x cot ) + y) , 

2hn 

cos(ft)' 

Solution (12.71) is defined over the finite time interval < t < 2-r j j . Note that particle 
trajectories have the following form (after time elimination) 

(x - A) 2 + (y- Bf = R 2 , 

where 

A = x /2 + u /f, B = y /2 + v /f, R 2 = (-x /2 + u /f) 2 + (y /2 - v /f) 2 . 

Here xo and yo are coordinates of the particle at t — ir/f. Thus, each fluid particle 
moves around a circumference. However, these trajectories do not close because solution 
becomes invalid as time t tends to or 2n/f. This solution describes essentially two- 
dimension fluid flow. In the case Mq + v q = it reduces to the rotationally symmetric 
solution in the polar frame of reference (13. 2p 



2 \ 2 r 2 ' 1 - cos(/t) 

Inversely, let us take solution u — 0, v — 0, h — ho of the RSW equations (jl.ip . 
corresponding to the state of rest. Transformation (12.61) yields the following rotationally- 
symmetric solution of the SW equations (12. 4p 

TT f^ 1 " T/ f 1 " U k ° 



1 + fH 2 ' 1 + fH 2 ' l + f 2 t 2 

(primes are omitted here). 

Using analogy with gas dynamic, we can conclude that these solutions belong to 
the class of "barochronous" [T9l 120] motions (as depth h depends only on time t). 



3. Finite transformations corresponding to the "complicate" generators 

It is well known that generator X = £, l d x i can be associated with one-parameter group of 
transformation G\ given by the finite relations of the form x l = ^(x 1 , x N , a), where 
a is a real-valued parameter. To determine these relations we shall solve the system of 
Lie equations 

— = C(x\...,x N ), x% =0 = x* (* = i,..., iV). (3.1) 

Let us find the finite transformations corresponding to the nontrivial infinitesimal 
generators Y 7 , Y 8 and Y" 9 . For the sake of convenience we shall use polar coordinates 
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(r,8), radial U and circular V velocity vector components related to (x, y) coordinates 
and (u, v) velocities, so that 

x = rcos8, y = rsin6, u = U cosd — V sin6, v = U sind + V cosd. (3-2) 

In these variables admissible symmetries Y 7 , Y s , Y 9 take the following form 

Y 7 = f~\l - cos(ft))d t + 2- 1 rsm(ft)d r - 2~\l - cos(ft))d e - 

-2- 1 (f/sin(/t) - frcos(ft))du - 2~\V + fr) sm(ft)d v - hsm(ft)d h , 

Y 8 = -Y 7 + 2f- 1 d t - d e , 

% = -2f- 1 sin(ft)d t - r cos(ft)d r + sm(ft)d e + 

+{Ucos{ft) + frsm{ft))du + (V + fr) cos{ft)d v + 2hcos(ft)d h . 

Integration of the Lie equations (13. 1ft associated with generator Yg gives the following 
result (t 7^ {2n + l)n/ f, n — integer): 



- 2 . , . , la(l + r 2 ' 

t = — arctan(ar) + r 



f v " y ' ^ ^ V l + a 2 r 2 

8 = 6 + arctan(r) — arctan(ar), 



fr (a 2 - 1)t \ / 1 + « 2 r 2 (3.3) 
' 2 l + a 2 r 2 ) \ a(l + T 2 Y 



l + a 2 r 2 J\ja(l + T 2 y a(l + r 2 )' 

Here r = tan(y), a = exp(— 2a), x(t) = 2irk/ f, where integer k is such that value t 
specified by initial condition at a = belongs to interval ((2k — 1)tv/ f, (2k + 1)tt/ f). In 
the case t = (2n + l)vr// solution of the equations (13. ip have the following form 

t = t =(2n + l)rc/f, f = r/V5, 6 = 6, 

U = Uy/^, V = (V+(a-l)(2a)- 1 fr)^, h = ah. 

Let us show that functions (13. 3p can be defined at points t = t* = (2n + 1)tt / f 
by continuity in compliance with ( 13. 4ft . The derived mapping will be continuously 
differentiable with respect to all its arguments. Let us investigate one-sided limits 
in formula (|3 .31) at t — ► t* =F (the remaining variables are fixed). In this case 
r(t) = tan(y) — > ±oo. Discontinuous at point t = i* function j arctan(ar(t)) has 
the following limits: rr/ f on the left-hand side and —n/f on the right-hand side of the 
discontinuity. Function x(t) is piecewise constant and possesses the value of 27m/ / on 
the left-hand side and 2w(n + 1) //on the right-hand side of point t = t*. Thus, function 
i(t, a) = j arctan(ar) + x(t) with t — > =F has one-sided limits which coincide and 
are equal to (2n + 1)tt/ f. This agrees with the first formula in (13 .4p . Correspondence 
of the remain functions in formulae (13.31) and (13.41) is easily established by calculation 
of their one-side limits at t = t*. Moreover, straightforward analysis shows that the 
mapping (13.31) redefined at points t = t* = (2n + l)vr// in compliance with formulae 
(13. 4p is continuously differentiable with respect to all the arguments. 
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- 2 / \ , v _ I r 2 + l 

t — — arctan(r + a) + X\t)i r 



6 = 6 + arctan(r) — arctanfr + a) 



u= {u+ fr(r 2 + ar-l)a\ /( r + a) 2 + l 



r+a 2 +l / V r 2 + l 



y= fr(2r + a)a_\ + (r + a)- + 1 



(3.5) 



2 (r + a) 2 + I J V r 2 + 1 ' r 2 + 1 

Here r and % are the same functions as before. 

Finite transformations corresponding to the generator Y 7 have the same form (13. 5p . 
where a = — cot(y) and Xi(t) stand for r and x(t)- Function Xi(t) is equal to 
(2k + l)7r//, where integer /c is such that value t specified by initial condition in the 
system (13.11) at a = belongs to interval (2kir/ f, 2(k + f). 

According to the general theory [5], Ej governing equations (11.11) presented in the 
polar coordinates (13. 2\i 

— + U— + —— - — -fV+— - 
dt dr r 86 r Or 

dV TT dV VdV UV £TT gdh n 

at or r 06 r r 06 

dh ld(rUh) ld(Vh) 

dt r dr r 86 
remain unchanged under the substitutions (13. 3p and (13. 5p corresponding to the finite 
transformations of the "complicate" generators. This fact allows one to formulate the 
following theorem based on the finite transformations (I3.3p . Similar theorem can be 
obtained using formulae (13.51) . 

Theorem 4 If a set of functions 

U = U(t,r,9), V = V(t,r,9), h = h(t,r } 6) 

satisfies the system of equations Ii3. 6]) then the following set of functions 

ff(a 2 -l)r 




U(t } r,6)=[U(i,f,6) + , 

2 cr + r z 

V(t,r,6) = (v(t, f ,S) - ^-^-°) \ . IJ*±<L. (3.7) 

h(t,r,6)= a * + /! 2 M ,f,S) 
a(l + t z ) 

also satisfies the same system. Here 

t = tan ( ^— ] , r = tan ( ^— ] , i = — arctan(ar) + 



2 J' f 



a(l + r 2 ) 

— — , 6 = 6 + arctan(r) — arctan(ar) 

1 + arr 2 



(3- 



Symmetries and exact solutions of the rotating shallow water equations 



10 



where a is an arbitrary positive constant, x(t) is the piecewise constant function which 
is equal to lukjf on the interval t G ((2k — l)n/f,(2k + l)7r//) ; k - - integer. We 
assume that t = t, f = r/y/a., 9 = 9 at points t = (2k + 1)71"//. 



4. Time-periodic exact solutions of the RSW equations 

The above-formulated theorem makes it possible to obtain new exact solutions of the 
RSW model by using known ones. Similar approach was applied to the gas dynamic 
equations in [211 122]. Obviously, equations (13. 6p have the following class of stationary 
rotationally-symmetric solutions 

U = U = 0, V = V(r), h = h(r) = ^- J (^— + fV^j dr + ho, (4.1) 

o 

where V(r) is an arbitrary smooth function, ho is a positive constant. According to the 
preceding theorem and formulae (13.71) the functions 

fr (a 2 - l)i 

i > 

fr (a - l)(ar 2 - 1) 
Y 1 + a 2 r 2 ' ( 4 - 2 ) 

= tan ^- 

is the solution of the equations (13. 6p . The principal prerequisite to the existence of 
periodic solutions of the form (14. 21) is the presence of nonzero Coriolis parameter /. Let 
us consider in detail two typical solutions from class (14. 2p . 

4-1. Example #1: Pulsation of the liquid cylinder 

Let us consider solution from class (14.11) corresponding to the state of rest: U = 0, 
V = 0, and h = h . By substitution of these functions into (14.2j) . the time-periodic 
exact solution of the RSW equations (13.61) 

fr (a 2 - l)r _fr (a - l)(ar 2 - 1) 

2 1 + a 2 r 2 ' " 2 1 + a 2 r 2 ' u ^ 

a(l + r 2 )h { " 




1 + a 2 r 2 



is obtained. Here, as before, r = tan(y). Particle trajectories are defined by the system 
of equations 

Tt =U, _ = -, rU = r , O\„ = 0o (4.4) 
and in this case have the following form 



<l i a 2 T 2(f) 

r(t) = r J — 9(t) = 9 + arctan(ar(t)) - arctan(r(t)). (4.5) 

l + T Z (t) 
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Formulae (j4.5p can be rewritten in a simpler form 
(x - A) 2 + (y - Bf = R 2 , 

where 

A = 2 _1 (a + l)r cos# , B = 2~ 1 (a + l)r sin 6 , R = 2~ 1 (a - l)r . 

Thus, each fluid particle moves round a circumference of radius R, which is proportional 
to the initial distance of the particle from the origin of the system of coordinates. 
All particles return to their initial positions at time T = 2ir/f, that is period of 
solution ( 14. 3D . The solution can be interpreted as pulsation of the liquid cylinder with 
impermeable boundary moving in accordance with the first formula (14.51) . where r is 
the initial cylinder boundary radius. The height of the liquid column, depending only on 
time, is given by last formula in (14.31) (ah is fluid surface height at initial time t — 0). 
We also note that time-periodic solution ( 14.3}) is the solution with constant potential 
vorticity Q = f/ho. 

4-2. Example #2: Pulsation of the liquid "drop" 
Let us consider the following solution from class ( 14.11) 

U = 0, V(r) = lr 2 , h(r) = £ ( r 4 + %r* + £Y (4.6) 



Ag \ 31 3Z 4 , 

where I = -f 2 ^a/(\2g) < 0. In this case the fluid is localized on the rotating plane 
within the cylinder < r < —f/l (we suppose, that / > 0). The depth of fluid h is 
maximal at r = and vanishes at r = —f/l. According to formulae (14.21) and (14.61) 
we arrive to the new time-periodic exact solution of the RSW equations ( 13. 6ft . which 
describes pulsation of liquid volume under the influence of gravitation and rotation. 

Integrating the equations of trajectories ( 14.41) with the functions U, V defined by 
(14, 2p and (I4.6p . we obtain the following result 



r(t) = r i 



l + a 2 r 2 (t) 



l + T 2 (t) ' (4.7) 

6(t) = (2C/ -1 + 1) arctan(ar(t)) - arctan(r(i)) + C x(t) + . 

Here C = lro^/a^, x(t) is a piecewise constant function which takes the value 2/nh/ f for 
t from the interval ((2k — 1)tt/ f, (2k + 1)tt/ f), where k is an integer. Constants ro and 
#o specify the position of the particle on the plane at initial time t = 0. 

Figure [1] presents the fluid depth depending on the radius r at different time 
moments t = 7m/(2f) with n — 0,1,2 (lines 0, 1, and 2, correspondingly). The fluid 
depth takes zero value at the circumference of radius 



/ y a(l + r 2 (t))' 

The graphs shown in figure [1] and others are plotted for a = 2, / = 1, g = 1. As the 
parameters change, the qualitative view of the figure remains the same. 
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0, 1, and 2, correspondingly). 



Formulae (14. 71) imply that all particles, which belong to the the circumference of 
radius r = f| (here M > m > are integers) at initial time t — 0, have 

closed trajectories. At £ = l-nMjf the particles return to their starting position at 
t = 0. Generally, trajectories of fluid particles are quasi-closed. This means that for any 
specified value e > it is possible to choose t £ > 0, so that the particle at t — t E will 
be located at a distance of no more than e from its initial position at t — 0. It should 
be pointed out that within the period of the solution T = 27r// only those particles, 
which belong to the boundary r = R*(t), where fluid depth h vanishes, return to their 
starting positions. In figure [2] closed trajectories of particles, located at t = in the 
circumference of radius r = l/(2v3) (curve 1) and r = l/y3 (curve 2), are shown. The 
graphics correspond to the choice (m,M) = (1,6) and (1,3), with the above-indicated 
values a, f and g. Arrows show the direction of the particles motion. 

Figure [3] shows the evolution of the material curves, which consist always of the 
same fluid particles. In contrast to the previous solution ( 14. 3ft . here an arbitrary closed 
material curve, except for circumferences with centre at origin of coordinates, never 
returns to its initial position. Propagation of the closed material curve (which is initially 
the circumference of radius r = 0.3 with its centre at point (x,y) = (0.4,0.5)) is shown 
in figure G^a) at the instants t = nn/(2f) with n = 0, ...,4 (curves 0-4). The trajectory 
of the fluid particle, initially located at point (x,y) = (0.7,0.5), is shown in figure [3]^a) 
by dashed line. As time growth this material curve transforms into helix. Figure [3^6) 
presents two material curves which form circumferences at t — (labeled as 1 and 2); 
and these material curves at the instant t — lir/f (helixes 3 and 4, correspondingly). 

Remark 1 The deformation of the material curve, shown in the Fig. is typical for 
all solutions from the class ffJfy where function V is not identically zero. Indeed, we 
can integrate equations of trajectories &4-4\ ) and represent the solution in the form fl^. 7| ) 
with C = (r y / «) _1 V(ro\/a). 
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Figure 3. The evolution of the material curves: a — short time behaviour; b — long 
time behaviour. 



5. Rotationally-symmetric flows 



The infinitesimal symmetries (I2.1j) -( l2~3l) allow us to construct invariant and partially 
invariant solutions to the RSW model. This section is devoted to a set of rotationally- 
symmetric submodels, which can be reduced to ODEs. These submodels are derived 
by invariant reduction of the RSW equations (13.61) with the use of two-dimensional 
parameterized classes of the optimal system of the subalgebra QL g . Therefore, the 
obtained solutions are essentially different, that is one solution does not reduce to the 
other by change of the variables. Note that these solutions can be obtained from the 
corresponding invariant solutions of the SW equations (12.41) using transformation ( 12.61) . 



(5-1) 



5.1. Stationary rotationally-symmetric flows in a ring 
Submodel 

UU r - r~ V 2 - fV + gh r = 0, 

U(V r + r~ x V + f) = 0, (rUh) r = 0. 

characterizes the class of stationary rotationally-symmetric solutions. Equations ( 15.11) 
are derived according to the invariant reduction of the RSW model (13.61) with the use of 
subalgebra (Y5, Y? + Ig) of admitted infinitesimal symmetries (in the polar coordinates 
the generators have the following form Y5 = d$ and Y? + Y$ — 2f~ 1 dt — de). 

Assuming U = 0, integration of (15. ip yields (14. ip . In case U 7^ 0, we obtain the 
following result 



U 



where C,- 



Ch v= C2_fr 
rh' r 2 ' 

1, 2, 3) are constants and 



F(h,r) = h 3 + pi(r)h 2 + V2 (r) = 0, (5.2) 



(pi(r) = - 
9 



1 ( Pr 2 



C? 
2T 2 



cl 

2gr 2 
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Figure 4. Stationary rotationally-symmetric flow in a ring bounded by characteristics: 
a — free surface height (line 1 corresponds to supercritical flow, line 2 — subcritical 
flow); b — velocity field (supercritical flow); c — velocity field (subcritical flow). 



The necessary condition of existence of solutions to the equation F(h, r) = with h > 
is inequality C\ > |C 2 |/2, which guarantees that (fx < for some values r. 

The straightforward analysis shows that equation F(h, r) = has two branches of 
solutions h = h(r) defined within the interval r* < r < r* (see figure Hta)). At the 
endpoints r = r* and r = r* the derivative vanishes: F^ = (3h + 2(fi)h = 0. Hence, 
function h{r) at the endpoints is finite h = h c (here h c = — §<^i(r*) or h c = — |^?i(r*)), 
whereas its derivative is unbounded: /i'( r *) = h'(r*) = —F r /Fh = oo. The relation 
U 2 /2 + gh = —gipiir) is fulfilled for the considered solution (15.21) . so that we get equality 
U 2 = gh in circumferences of radii r = r* and r = r*. Thus, the surfaces r = r* and 
r = r* are sonic characteristics of the model. Recall that the surface $(£, r, 9) = const is 
the sonic characteristic of equations (I3.6P on the solution (U,V,h) if function $(t,r, 6) 
satisfies to equation $ t + U& r + r~ x V^0 = ±^fgh. 

The flow is supercritical (subcritical) in case Fr > 1 (Fr < 1), correspondingly. 
Here Fr = q/y/gh is Froude number and q = \J\J 2 + V 2 is the magnitude of the 
velocity vector. It follows from ( 15. 21) that q 2 = 1C\ — C 2 / — 2gh. Hence, Fr = 1 
at h = h s = (2Ci — C*2/)/(3(7). One can easily prove that h c < h s . Therefore, the lower 
branch of solution to equation F(h, r) = corresponds to the supercritical flow. In 
figure H](a) are presented the lower (line 1) and the upper (line 2) branches of solution 
to equation F(h, r) = 0. In figures H](6) and H](c) the velocity fields (w, v ) are shown 
for lower branch of free surface height h = h(r) and upper branch, correspondingly. 
The graphics are derived for Cj = 1, g — 1, / = 0.1. As the parameters change, the 
qualitative view of the graphics remains the same. 

It is interesting to note that the class of stationary rotationally-symmetric solutions 
of the RS W equations corresponds to the class of solutions of the SW model which 
is invariant under rotation (generator Z^) and nontrivial projective transformation 
(generator Z 7 + Z 8 ). On the basis of projective transformation exact invariant solutions 
of the S W model and two-dimension gas dynamic equations (with special state equation) 
were obtain and studied in [21], [2H [23] and others. 
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5.2. Collapse of a liquid ring. Regime #1 

Let us consider the submodel derived by the invariant reduction of the RSW equations 
( 13. 6ft using two-dimensional subalgebra (Y&, Yj) of infinitesimal transformations. 
According to the algorithm of invariant solutions construction, let us solve the equations 

Y 5 J = 0, Y 7 J = 

to obtain the set of basic invariants 

f l-cos(Jt) fr* Mft) fr* , 

J ~\ r 2 > rU 2 l-cos(/t)' + 2 ' 
The representation of an invariant solution is 

u=m + fi cot (fii v= m* h= m, (5 . 3) 

r 2 V 2 / r 2 r z 

where A = (1 — cos(/t))r~ 2 . By substituting the representation of solution into (13. 6p . 
one obtains the system of ODEs for unknowns (p, if), rj 

(\(<p 2 + 2gr ] ))' + if) 2 = 0, A^' = 0, (A^)' = 0. (5.4) 

By virtue of the second equation ( 15.4D either ip = or if) = const. In the first case 
we get the following class of solutions 

cp = o, if; = if>(\), V = j (\oVo ~Y g S x ^ 2(Z/) dU ) • (5 ' 5) 

Here is an arbitrary smooth function, and rjo, Ao are constants. Solution ( 15. 3ft . 
(15. 5p is defined over the finite time interval t e [to,27r//), t > 0. It has a singularity 
at r — > oo. In some cases, depending on the choice of function ■0, the solution has 
singularity at r — > 0. 

Notice, that surfaces A(t,r) = const are contact characteristics on the solution 
( 15. 3p . (15.51) . This allows one to interpret the solution as a collapse of a liquid ring (or a 
liquid cylinder) R\(t) < r < compressed by pistons which move according to the 

law 



Ri(t) = R i0 sin(/t/2) (0 < R w < R 



20 ; 



so that the pistons surfaces coincide with the contact characteristics. In the course of 
the collapse of a liquid ring (with t — > 2n/f), the radial velocity U = v cot (y) in 
the range of Ri(t) < r < R2(t) is bounded, whereas the liquid depth h as well as the 
circular velocity V (with any ip ^ const) increases infinitely. 

The other class of solutions of submodel ( 15. ip is derived in case ip ^ 0, if) = C2 = 
const. Whereas, 77(A) = C^/{\(p) and function y?(A) are defined from cubic equation 
V? 3 + (Cf - Ci/X)(p + 2gC 2 /X = 0. The solution of the RSW equations (EEJ is obtained 
by using representation of solution (15. 3p . 

These unsteady solutions of the RSW equations can be also derived from well-known 
stationary rotational-symmetric solutions of the SW equations using transformations 
(12. 6p which relate the models. 
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5.3. Collapse of a liquid ring. Regime #2 
Submodel 

V = + - V? 2 - 2077, ^' = -(2^ + r]' = -Aiprj. (5.6) 

arises from the invariant reduction of the RSW equations (13.61) with the use of the 
subalgebra (Y 5 ,Y e ). Indeed, in polar coordinates (13.21) generator F 6 has the following 
form: Y 6 = rd r + Udu + Vdy + 2hdh- The set of basic invariants of the generators Y 5 
and Yq is 

J= {t,r- x U,r- x V,r- 2 h}. 

Thus, we get the following representation of the invariant solution 

U = r<p(t), V = ii)(t), h = r 2 r](t). (5.7) 

Substitution of representation (15 .7p into the RSW model (13.61) produces ODEs (15. 6p . 
In the general case it is difficult to get a solution of the submodel (I5.6P in closed form. 
Let us take ip(t) — — //2. Then the second equation in (15.61) is fulfilled automatically, 
whereas the solution to the first and the third equations is given in an implicit form 

v = m = ±^ m - £ + (»* - 2ff » + £) y/i, *w = 4£4^. (5 ' 8) 

Here arbitrary constants <^ and T]o > are values of functions tp(t) and r)(t) at t = 
(initial data for ODEs ( 15. 6ft ). Function (p(t) is strictly decreasing by virtue of the first 
equation (15.61) (due to ip — —f/2). It follows from the last equation (15.61) that function 
rj(t) increases as time grows if <p(t) < 0, otherwise rj(t) decreases. 

If constant ipo < 0, the minus sign precedes the root in the first formula ( 15. 81) . 
In fact, (p(0) = ip(r] ) < and (p(t) is strictly decreasing function. It is evident from 
the second formula (15.8P that t^T*>0as?7^oo (integral in (15.81) converges). 
The solution (15. 7p . (15.81) can be interpreted as collapse of a liquid ring during the 
finite time interval [0, T*]. At t = liquid is located in the ring (or in the cylinder) 
< -Rio < r < -^20- As time grows, impermeable boundaries of the ring move in 
accordance with the law 

r = R t (t) = Riir}) = R i0 (vo/v) 1/4 , 
providing the fulfillment of boundary condition R'^t) = U(t,Ri(t)). In contrast to the 
previous solution ( 15.31) . (15.51) at t — » T* the circular velocity V vanishes, whereas the 
radial velocity U and the depth h infinitely increase. At that, the moving pistons collapse 
into the origin r = 0. 

Let us take <po > 0. In view of <^o — v(0) = 0(f]o) > 0, the plus sign precedes 
the root in first formula ( 15. 8H . As mentioned above, function rj decreases while p > 0. 
It is obvious from (I5.8P that radicand vanishes at some value rj = rji < tjq. At time 
t\ — t(j]i) < oo function cp(t) changes sign from positive to negative and function rj(t) 
grows for t > t\. The behaviour of solution for t > t% is similar to the previous case 
fo < 0. Thus, the solution can be interpreted as a partial spreading (0 < t < ti) and 
collapse (ti < t < T«) of a liquid ring. 



Symmetries and exact solutions of the rotating shallow water equations 



17 



6. Conclusions 

In this paper we have studied symmetry properties and some classes of exact solutions of 
the RSW model using group analysis. We have determined the 9-dimensional Lie algebra 
of admissible infinitesimal generators. We have shown that the derived Lie algebra of 
symmetries is isomorphic to the Lie algebra of infinitesimal transformations admitted 
by (2+l)-dimensional SW equations. This allows one to employ its known optimal 
system of subalgebras for construction of essentially different invariant solutions of the 
RSW equations. Moreover, we have found change of variables which reduce the RSW 
equations to the SW model. This transformation allow one to construct and study 
solutions of the RSW equations using solutions of the SW model and vice versa. We 
have also derived and analysed finite transformations corresponding to the nontrivial 
symmetries of the RSW model. Using these transformations we have generated new 
time-periodic exact solutions of the RSW equations. These solutions describe fluid 
flow with quasi-closed particle trajectories and may be interpreted as pulsation of 
liquid volume under the influence of gravity and Coriolis forces. Based on some two- 
dimensional parameterized classes of the optimal system of subalgebras we have reduced 
the RSW model to the ODEs and integrate them. In particular, we have constructed and 
studied exact solutions of the RSW model describing rotational-symmetric stationary 
flows in a ring bounded by characteristics as well as time dependent flows which describe 
spreading and collapse of a liquid ring. We have also pointed out exact solutions of the 
SW equations corresponding to the obtained solutions of the RSW equations. 
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